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ABSTRACT 



Or 



Galaxy dynamical masses are important physical quantities to constrain galaxy evolutionary models, especially at high redshifts. 
However, at z > 2 the limited signal to noise ratio and spatial resolution of the data usually do not allow spatially resolved kinematical 
modeling and very often only virial masses can be estimated from line widths. But even such estimates require a good knowledge of 
galaxy size, which may be smaller than the spatial resolution. Spectroastrometry is a technique which combines spatial and spectral 
resolution to probe spatial scales significantly smaller than the spatial resolution of the observations. Here we apply it to the case of 
high-z galaxies and present a method based on spectroastrometry to estimate dynamical masses of high z galaxies, which overcomes 
the problem of size determination with poor spatial resolution. We construct and calibrate a "spectroastrometric" virial mass estimator, 
modifying the "classical" virial mass formula. We apply our method to the [O III] or Ha emission line detected in z ~ 2 - 3 galaxies 
from AMAZE, LSD and SINS samples and we compare the spectroastrometric estimator with dynamical mass values resulting from 
full spatially resolved kinematical modeling. The spectroastrometric estimator is found to be a good approximation of dynamical 
masses, presenting a linear relation with a residual dispersion of only 0.15 dex. This is a big improvement compared to the "classical" 
virial mass estimator which has a non linear relation and much larger dispersion ( 0.47 dex ) compared to dynamical masses. By 
applying our calibrated estimator to 16 galaxies from the AMAZE and LSD samples, we obtain masses in the ~ 10 7 - 10 10 Mq range 
extending the mass range attainable with dynamical modeling. 



oo 
o 



- 1—^ 

X 



Key words. Methods: data analysis - Techniques: high angular resolution - Techniques: spectroscopic - Galaxies: fundamental pa- 
rameters - Galaxies: high-redshift - Galaxies: kinematics and dynamics 

1. Introduction but there remain many dark areas because of the lack of 

observations that can constrain the complex p hysical processes 
The dynamical properties of galaxies have a fundamental role involved. In galaxy evolutionary models (iBlumenthal et all 
^ . in the context of galaxy formation and evolution. They are a 19841 iDavis etall fl985t ISmingel et alJ 120061 iMo et al.1 119981 
0j ' key prediction of theoretical models and are the most direct way mer gers are believed to play an important role for galaxy 
to probe the content of dark matter. In particular, measuring f ormat ion and evolution but there is also observational evidence 
the dynamical mass of a galaxy is the most direct way to for the existence o f rotating disks with high s t ar formation 
constrain the mass of the dark matter haloes, and this quantity rates at z ~ 2 (Forster Schr eiber et al.l 2009, C resci et al.l 
is a fundamental prediction of cosmological cold dark matter 12009), Whi i e there are many dynamical studies on extended 
models of hierarchical structure formation. In recent years, our samples of z ~ 2 object s dGenzel et al.ll2006t Idenzel et al.l|200~ 
observational knowledge in this field has increased enormously, poster Schreibe retal] l2006bt iForster Schreiber et alJ 1200 
ICresci et all 12009; Erb et al.l 120061) . little is known about the 



Send offprint requests to: gnerucci@arcetri.astro.it dynamics of galaxies at redshift z > 2.5 and beyond where 



Based on observations collected with European Southern only a handful of objects have been observe d (Nesvadba et al 

Observatory/Very Large Telescope (ESO/VLT) (proposals 075.A- | 2006|; iNesvadba et all 120071 : [Nesvadba et a ll 120081: iJones et aL 

0300, 076.A-071 1 and 178.B-0838), with the Italian TNG, operated by 201Q : Law et all 120091: iLemoine- Busseroll e et all 12009: 

FGG (INAF) at the Spanish Observatorio del Roque de los Muchachos, ISwinbank et alJ 120071: ISwinbank et alJ 120091) . The importance 

and with the Spitzer Space Telescope, operated by JPL (Caltech) under f studying galaxies at very high redshifts (z ~ 3 - 4) and 
a contract with NASA. 
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measuring their dynamical masses can be summarised as 
follows: this redshift range is before the peak of the cosmic 
star f ormation rate (see, for example, iHopkins & Beacoml 
2006: IMannucci et all l2007h . only a small fraction (~ 15%, 
Pozzetti et alj ' 2007F of the stellar mass has been created, it is 
in this redshift range that the most ma ssive early-type gal axies 
are expected to form (see, for example, Sar acco et al.ll2003h and 
the number of galaxy mergers is much larg er than at later times 
dConselice et all 120071: IStewart et all 12008). As a consequence 
the predictions of different models start to diverge significantly 
at z > 2.5, and this divergence can be very large at z ~ 3 - 4. 

We have recently measu red dynamical masse s in a sizable 
samp le of galaxies at z ~ 3 dGnerucci et al.ll20lll , ICresci etldl 
12010) . One of the lesson learned from our own work and from 
the literature is that the data for objects at such high redshifts 
often suffer from poor signal to noise (hereafter S/N) which does 
not allow spatially resolved kinematical studies and complete 
dynamical modeling. For this reason one can only estimate the 
dynamical mass of a galaxy by applying the virial theorem to its 
integrated spectrum. 

Virial mass estimates are affected by large systematic errors. 
Apart from those due to the unverified assumption that the sys- 
tem is virialized, one of the principal problems is the estimate of 
the size of the galaxy, which often suffers from the low intrin- 
sic spatial resolution. In this paper we present an alternative to 
the classical virial mass estimate based on the technique of spec- 
troastrometry with the aim of obtaining a more a ccurate dynam- 
ical m ass estimator. We follow up on the work of iGnerucci et al.l 
(2010) about the application of the spectroastrometry technique 
to the study of the dynamics of rotating gas disks. 

This paper is based on the data pertaining to two projects 
focussed on studying metallicity and dynamics of high- 
redshift galaxies: AMAZE (Assessing the Mass- Abundance red- 
shift Evolution) dMaiolino et al.ll2008al iMaiolino et al] 12008b) 
and LSD (Lyman- break galaxies Stellar populations an d 
Dynamics) (IMannucci & Maiolinol2008l IMannucci et al1l2009l) . 
Both projects use integral field spectroscopy of samples of z ~ 
3-4 galaxies in order to derive their chemical and dynami- 
cal properties. In both projects we make use of data obtained 
with the Spectrograph for Integral Field Observations in the 
Near Infrared (SINFONI) at the Very Large Telescope (VLT) of 
the European Southern Observatory (ESO). Integral field spec- 
troscopy has proven to be a powerful tool to study galaxy dy- 
namics as it provides two dimensional velocity maps, without 
the restrictions of longslit studies, which are also plagued by un- 
avoidable light losses and dynamical biases due to a possible 
misalignment between the slit and the major axis. 

The AMAZE sample consists of ~ 30 galaxies selected with 
deep Spitzer/IRAC photometry (3.6 - 8/im), an important piece 
of information to derive a reliable stellar mass. These galaxies 
were observed with SINFONI in seeing-limited mode. The LSD 
sample is an unbiased, albeit small, sample of LBGs with avail- 
able Spitzer and HST imaging. SINFONI observations were per- 
formed with the aid of adaptive optics in order to improve spatial 
resolution since the aim of this project was to obtain spatially re- 
solved spectra for measuring kinematics and gradients in emis- 
sion li nes. 

In IGnerucci et al.l d201 ll) we performed a kinematical anal- 
ysis of 33 z ~ 3 galaxies from the AMAZE and LSD sam- 
ples. We found that ~ 30% of the objects show gas kinemat- 
ics consistent with that of a rotating disk. For these "rotat- 
ing" galaxies we performed kinematical modeling by fitting ro- 
tating disk models and found dynamical masses in the range 
~ 1.8 x 10 9 M© - 2.2 x 10 n M©; we also found that for the ma- 



jority of the objects the contribution of turbulent motions is im- 
portant compared to the ordered motions, suggesting that most 
rotating galaxies at z ~ 3 are dynamically hot disks. We then 
built the baryonic Tully-Fisher relation at z ~ 3. Our measure- 
ments indicate that z ~ 3 galaxies have a stellar-to-dynamical 
mass ratio smaller than the value in the local universe, confirm- 
ing the redshift evolution of the relation already found at z ~ 2. 
However, the large scatter of the points may also suggest that, at 
this redshift, the relation is not yet in place, probably due to the 
young age of the galaxies. 

In Sect. [2] we briefly summarize observations and data re- 
duction for the data used in this paper. In Sect. |4] we describe 
the principle of the spectroastrometry technique and its ap- 
plication to local blac k hole mass measurements presented in 
IGnerucci et al ] (120101) . In Sect. we introduce the "classical" 
virial mass estimator and calibrate it by comparing it with the 
more accurate dynamical mass estimates. In Sect. |5]we present 
our spectroastrometric dynamical mass estimator and in Sect. 15. II 
we calibrate it by comparing it with the more accurate dynamic 
mass measurements. In Sect. [6] we apply our method to galaxies 
in the AMAZE and LSD samples and discuss the results. Finally, 
in Sect.|7]we draw our conclusions. In Appendix lAl we use sim- 
ulations to understand how our estimator is affected by various 
dynamical or instrumental parameters. 

In the rest of the paper, we will adopt a KCDM cosmology 
with H = 10kms-\ D. m = 0.3 and Q A = 0.7. 



2. The data 

Complete descriptions of the AMAZE and LSD programs, of 
their sam ple selection, ob s ervatio ns an d data reduction are pre - 
sented in Mai olino et"ai1 (l2008al) and IMannucci et al.l J2009). 
Here we present a brief summary of the observations and data 
reduction. 

The near-IR spectroscopic observations were obtained 
using SINFONI, the integral field spectrometer at VLT 
(lEisenhauer et al.ll2003l) . For AMAZE galaxies, SINFONI was 
used in its seeing-limited mode, with the 0.125" x 0.25" pixel 
scale and H+K grism, yielding a spectral resolution R ~ 1500 
over the spectral range 1.45 - 2.41/im. For LSD galaxies, 
SINFONI was used with the Adaptive Optics module using a 
bright star close to the galaxy to guide the wavefront correction 
system. The (K-band) seeing during the observations was gener- 
ally better than 0.8". Data were reduced with the ESO-SINFONI 
pipeline. The pipeline subtracts the sky from the temporally con- 
tiguous frames, flat-fields the images, spectrally calibrates each 
individual slice and then reconstructs the cube. Individual cubes 
were aligned in the spatial direction using the offsets of the po- 
sition of the [O III] or Ha line emission peak. The atmospheric 
absorption and instrumental response were taken into account 
and corrected by dividing with a suitable standard star. The gas 
kinematics used in the following analysis were obtained from 
[O III] AA 4959, 5007 A and Ha for galaxies at z ~ 3 and z ~ 2, 
respectively. 

2. 1 . Complete dynamical modeling of the objects 

In lGnerucci et al.l d201 lb we presented the spatially resolved dy- 
namical modeling of the AMAZE and LSD objects. In this paper 
we will use the results of that modeling as a basis to calibrate the 
dynamical mass estimator we are developing. We refer to that 
paper for any detail of the dynamical modeling procedure. 
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From the AMAZE and LSD samples we selected a subsam- 
ple of galaxies that show a velocity field consistent with a ro- 
tating disk. In iGnerucci et al.l 1201 ll we implemented a simple 
but quantitative method to identify smooth velocity gradients, 
which are the signature of a gas kinematics consistent with disk 
rotation. This method is based on classifying a galaxy as "ro- 
tating" or "non rotating" according to whether its velocity map 
is well fitted (or not) by a plane. Using this method we selected 
from the AMAZE and LSD samples 1 1 "rotating" objects for the 
complete dynamical modeling, but we also identified 6 objects as 
"unclassified" because of the low spatial res olution or S/N (see 
section 3 and table 2 of Gnerucci et al. I l201lh . For the "rotating" 
objects we performed a complete dyna mical modeling. 

The adopted dynamical model (see IGnerucci et al]|201 ll for 
details) assumes that the ionized gas is circularly rotating in a 
thin disk neglecting all hydrodynamical effects, thus the disk 
motion is entirely determined by the gravitational potential. The 
galaxy gravitational potential is generated by the galaxy mass 
distribution that is modeled by an exponential disk mass distri- 
bution. The velocity along the line of sight for a given position 
on the sky is then derived by taking into account geometrical 
projection effects. The model takes also into account instrumen- 
tal beam smearing and binning over detector spatial and velocity 
axes. 

The spatial resolution of the final data cubes was estimated 
follo wing kinematical co nsiderations and was typically ~ 0.6" 
(see lGnerucci et all201 ll for details). 

The values of the model parameters (e.g. position angle of 
the disk line of nodes, inclination of the disk, characteristic ra- 
dius of the disk, dynamical mass) are obtained by fitting this 
model to the object velocity map. In particular it yields the best 
fit value for the dynamical mass of the galaxy Md v „, defined 
as the total mass enclosed in a 10 kpc radius (e.g. [Cresci et"al] 
2009) and this is the dynamical mass value we will compare 
to our estimator prediction. In table Q] we report the dynamical 
masses of the best fi t models for the various objects analyzed in 
IGnerucci et ail d201 ll) . 

3. Virial mass estimator 

The main and more reliable technique to estimate the dynami- 
cal mass of a high redshift galaxy is the modeling of spatially 
resolved gas kinematics. This technique requires that an object 
is well spatially resolved and has a high S/N to constrain all the 
model parameters with sufficient accuracy. Unfortunately these 
two requirements are not often met for high redshift galaxies, 
which are in many cases only marginally spatially resolved and 
with poor S/N. In these cases a method to estimate the dynam- 
ical mass of the object is to use the integrated spectrum of the 
object and apply the virial theorem. 

Applying the virial theorem, the dynamical mass enclosed in 
a sphere of radius R c is given by the simple following equation. 

AV 2 R c 

M(R C ) = f—Q± (1) 

where AV is the velocity dispersion along the line of sight which 
is usually estimated from, e.g., FWHM or velocity dispersion of 
the line profile or from the observed velocity gradient. The /- 
factor takes into account the geometry of the mass distribution 
and velocity projection along the line of sight; its actual value 
also depends on the choice of the velocity estimator. 

The choice of the appropriate / is based on estimating its 
average value over several galaxy mass distribution models, but 



also depends on how the gas velocit y and the galaxy characteris - 
tic radius ar e estim ated. Following Bi nnev & Tremairiel (12008). 
lEpinat et alj d2009t) use a value of 2.25, an average value of 
known galactic mass distributio n models. Taking into account an 
average inclination correction, lErb et al.l (2006) obtain a value 
of 3.4. Finally, in the case of dispersion-dominated objects, 
Forster Schrei ber et alJ {2009) adopt an isotropic virial estimate 
with a value of 6.7, appropriate for a variety of galaxy mass dis- 
tributionfl 

The principal problem of the virial mass estimate lies in the 
estimate of the characteristic dimension R c . First one needs a 
continuum image of the source with large enough S/N to esti- 
mate, e.g., the half light radius. This size has then to be cor- 
rected to take into account the finite spatial resolution of the ob- 
servations; indeed one is often dealing with objects which are 
marginally spatially resolved and the estimate of R c can lead to 
large uncertainties and systematic errors. Second, one needs to 
estimate the velocity dispersion of the gas from the integrated 
spectrum. By combining gas velocity dispersion with contin- 
uum sizes one has to assume that gas emission is co-spatial with 
continuum emission and, in general, this assumption cannot be 
straightforwardly verified. Indeed, there is some disagreement 
in the few high-z studies regarding th e comparison of c ontin- 
uum and emission line morphologies. iKriek et all (2009) com- 
pared HST/NIC2 and SINFONI Ha morphologies of a spectro- 
scopic sample of massive galaxies at z ~ 2.3, and found in many 
cases significantly diff e rent morphologies. On the other hand 
Forster Schrei ber et alJ d201 ll) observed six galaxies at z ~ 2 
with HST/NIC2, finding a good agreement with SINFONI Ha 
mo rphologies. 

Forster Schreiber et a D (l2006bL (2009) introduce a variation 
to the classical virial mass estimator by using the projected cir- 
cular velocity at the turnover in eq.Q] They estimate this quantity 
by measuring the amplitude of the observed velocity gradient or 
the FWHM of the line in the integrated source spectrum. Then 
they multiply these quantities for an appropriate conversion fac- 
tor obtained by measuring their average ratio with the projected 
circular velocity at the turnover i n a subsample of galaxies mod- 
eled with rotating disks (see also lTacconi et al. 2006, 2008). 

Using the data from the AMAZE and LSD projects (see 
Sect.|2]for more detailed description) we can calculate the "clas- 
sical" virial mass for several galaxies (using eq.Q]) and compare 
it with the more accurate estimates f rom dynamic al modeling of 
spatially resolved gas kinematics in IGnerucci et al.l d201 ll) (see 
Sect. |23]). 

In particular we estimate the characteristic radius R c from the 
Gaussian HWHM of the emission line image of the galaxy (not 
having detected any continuum for most ob jects) correcting th is 
value for the instrumental beam smearing dBouche et al.l l2007). 

In the left panel of Fig. Q] we compare virial products (/j, V j r = 
FWHM 2 RJG) with masses resulting from full dynamical mod- 
eling (Mdyn) in IGnerucci et al.l d201 lb . The numbers plotted in 
the figure are also reported in table Q] We also extend this com- 
parison to a su bsample of 8 rotating disks from the SINS survey 
at z ~ 2 (F orster Schreiber et al. 2006 al iForster Schreiber et all 
l2009LlGenzeletalJl2008LICresci et alJl2009l) selected for their 
high S/N, "rotating disk" gas kinematics and robust dynamical 
modeling. 

For the comparison we do not use the dynamical mass Md yn 
itself, but (Md y „sin 2 i), where i is the inclination of the galaxy 
disk. Due to the coupling of dynamical mass and inclination, 

1 The / values presented here correspond to using the velocity dis- 
persion of the line as an estimate of velocity in eq. [T] 
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Fig. 1. Calibration of the dynamical mass estimator for the AMAZE, LSD and SINS objects. Left panel: calibration of the "classical" 
virial product yU,.,, = FWHM 2 R C /G. Dashed line: best-fit linear relation. Dotted lines: rms dispersion of the residuals from the 
fitted linear relation (0.3 dex). Dot-dashed line: best-fit non linear relation. See Sect. [3] for details. Right panel: calibration of the 
"Spectroastrometric" virial product fi spec = FWHM 2 r spec /G. Dashed line: best-fit linear relation. Dotted lines: rms dispersion of 
the residuals from the fitted relation (0.15dex). See Sect. l5.1l for details. 



Mdy„sin 2 i is the quantity which is less affected by systematic 
errors due to the uncertain disk inclination which is difficult to 
const rain especially with low S/N and barely resolved sources 
(see lGnerucci et all d201 lb ). Moreover, the uncertain disk incli- 
nation also affects the virial mass estimate because the FWHM 
is determined by the line of sight velocity dispersion. 

We note that the virial product yL vir and Md v „sin 2 i have a non- 
linear relation with large dispersion which is not the best cir- 
cumstance for a reliable prediction of Md m sin 2 i from yU l7 >. This 
is to be expected, to some extent, since /i,,,> depends on the line 
FWHM and Md yn depends on the velocity gradient. Therefore 
source size and seeing are both going to affect /z Vi > and Md yn in 
opposite ways : for a given velocity field, a reduction in spatial 
resolution will increase the spatially unresolved rotation and the 
line FWHM will get larger while, at the same time, the velocity 
gradient will be reduced. This might explain the measured slope 
in Fig.QJ. 

By imposing the expected linear relation log(jU V! >) = 
\og( Md yn sin 2 i) - log/', the mean rms dispersion of the fit resid- 
uals is 0.3 dex. By letting the slope vary freely, one obtains 
log(ju v! >) = (0.66 + 0.05)log(M £ / y „s/n 2 - log/' and the disper- 
sion reduces to 0.21 dex. The "classical" virial mass estimate 
can be biased by systematic errors principally originating from 
dimension estimates and therefore it is not a robust proxy for the 
dynamical mass. 

4. The Spectroastrometric technique 

Spectroastrometry, orig inally introduced by iBeckersi dl982h . 
IChristv et aD dl983h and lAime etall dl988h to detect unresolved 
binaries, has been used by several authors to study pre main se- 
quen ce binaries and the gas disks surrounding pre main sequence 
stars ([Baines et alJl2004t iPorter et~aT]l2004l 12005b iTakami et alJ 
120031 Iwhelan et alj2005l) . More recently. iGnerucci et alJ d2010) 
studied the application of spectroastrometry to constrain the 



kinematics of rotating gas disks in local galactic nuclei and to 
measure the mass of the putative supermassive black holes using 
either longslit spectra or integral field spectra. 

The fundamental advantage of the spectroastrometric 
method is that, in principle, it can provide position measure- 
ments on scales smaller than the spatial resolution of the obser- 
vations. This is due to the ability to separate the various spectral 
features of the source and obse rve the spatial l i ght d istribution 
of these features separately (see lGnerucci et af] d2010t) for a de- 
tailed analysis and discussion). 

4. 1 . Measuring masses with spectroastrometry 

Here we explain the basics of spectroastrometry and its ap- 
plication to measuring the masses of local Black Holes (here- 
after BH), i.e. the application to dynamical studies of gas disks 
rotating around poin t-like mass distributions, as described in 
IGnerucci etall d2010t) . 

The spectroastrometrical method consists of measuring the 
photocenter of light emission in different wavelength or veloc- 
ity channels. On a longslit spectrum of a continuum subtracted 
emission line (the so-called position-velocity diagram) the ro- 
tation curve denotes the mean gas velocity as a function of the 
position along the slit whereas the spectroastrometric curve pro- 
vides the mean positio n of the emitt i ng ga s as a function of 
velocity (see Fig. 1 of IGnerucci et al.l d2010l) ). The two curves 
clearly analyze the same spectrum from complementary points 
of view. 

IGnerucci et alJ (120101) performed an extensive set of simula- 
tions in order to understand how the spectroastrometric curve is 
affected by the object's own dynamical features or by the instru- 
mental setup. Such simulations show that the presence of a su- 
permassive black hole is revealed by a turn-over in the spectroas- 
trometric curve, with the high velocity component approaching 
a null spatial offset from the location of the galaxy nucleus. All 
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the relevant information about the BH is encoded in the high 
velocity range of the spectroastrometric curve, which is almost 
independent of the spatial resolution of the observations and of 
the intrinsic line flux spatial distribution. According to the sim- 
ulations, the use of spectroastrometry can allow the detection of 
BHs whose apparent size of the sphere of influence is as small 
as ~ 1 / 10 of the spatial re solution. 

iGnerucci et al.l ((2010) then provided a simple method to es- 
timate BH masses from spectroastrometric curves. This method 
consists of the determination of the spectroastrometric map, that 
is the positions, in each velocity bin, of the emission line pho- 
tocenters on the plane of the sky. These are obtained by com- 
bining longslit spectra centered on a galaxy nucleus at different 
position angles, but can be trivially obtained from integral field 
spectra. From the high velocity points in the spectroastrometric 
map one can then obtain a rotation curve to easily estimate the 
BH mass. This method has been tested with simulated data. The 
test demonstrates the possibility to reconstruct the rotation curve 
down to radii of ~ 1 / 10 of the spatial resolution of the data and, 
with seeing limited observations (~ 0.5"), to being able to de- 
tect a BH with mass 10 5 5 Mq(D/ 3.5 Mpc), where D is the galaxy 
distance, a factor ~ 10 better than can be done with the classical 
method based on rotation curves. 

The ~ 1/10 gain with respect to the spatial resolution clearly 
depends on the S/N of the data which ultimately determines the 
accuracy with which we can estimate photocenters. 

In this paper we make use of the results of IGnerucci et alj 
(1201 Oh to apply the spectroastrometric technique to the study 
of the gas dynamics in high redshift galaxies. In particular we 
use spectroastrometry to improve the "classical" virial mass es- 
timates. As noted in the introduction we work on IFU data 
(SINFONI integral field spectra). 

The principal difference with respect to its application to lo- 
cal BHs is that we deal with rotating gas disk dynamics driven by 
extended, instead of point-like, mass distributions. But in both 
cases spatial resolution is an issue, since the rotating disks are 
often barely resolved. 

5. Spectroastrometric mass estimator 

The spectroastrometric dynamical mass estimator presented in 
this paper is based on measuring the characteristic dimension 
of the object by means of the spectroastrometric technique. As 
for the classical virial mass estimates we assume that the ob- 
ject is a rotating disk. Following this assumption, we expect 
that the redshifted gas is located principally on one side of the 
disk, whereas the blueshifted gas is located on the other side. 
Therefore, if we ideally obtain images of the object in the "red" 
and "blue sides", these two images have to be spatially shifted 
due to the rotation. We estimate the characteristic radius of the 
object by measuring this shift between the red and blue sides 
with spectroastrometry. In the rest of the paper we will refer to 
"sizes" based on spectroastrometry. This is indeed a new mea- 
sure of size slightly different from all other measures typically 
used (e.g. half-light radius, scale length, HWHM, etc.). This 
"size" is in fact the average distance between the redshifted and 
blueshifted gas weighted on the surface brightness of line emis- 
sion. By means of simulations (see Appendix IA1 we can assess 
that for an exponential disk model this "spectroastrometric size" 
correspond on average to a fraction of ~ 0.75 of the disk scale 
length, although this relation depends on the assumed disk pro- 
file. 

The first step is to obtain a spatially integrated spectrum of 
the source. We then fit the observed line profile with a simple 



Gaussian function and estimate the line width (<x or FWHM). 
Using the fitted line profile model we select the "red" and the 
"blue" sides of the emission line. The red side is identified by 
the wavelength range bordered by the central wavelength of the 
line and the most extreme "red" wavelength for which all the 
following conditions are satisfied: (i) the line model is greater or 
equal to 5% of the maximum amplitude, (ii) the data spectrum 
is greater than or equal to the rms of the fit residuals and (iii) 
the wavelength is less than 3<x from the central wavelength. We 
similarly define the wavelength range for the "blue" side of the 
line (see Fig. [2j. We then collapse all the velocity planes in the 
red and blue sides of the line to obtain the "red" and "blue" im- 
ages of the object, respectively. The spatial plane related to the 
central wavelength bin is added to the red and blue images but 
with a weight given by the fraction of the bin lying on the red 
and blue sides (see Fig. |2j. 

Total spectrum of the object 
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Fig. 2. Example of the integrated spectrum for an object 
(SSA22a-C16). Overplotted is the fitted model of the line (red 
continuous line), the central wavelength (red dashed vertical 
line), the rms of residuals (two dotted horizontal lines) and the 
selected bins for the blue and red sides (respectively blue and red 
filled bins) and the central bin (brown filled bin). 

The next step is to fit such images with a two-dimensional 
Gaussian function and determine the position of the "blue" and 
"red" light centroids (see Fig. |4|. We can now define the spec- 
troastrometric characteristic radius (r spec ) as the half of the dis- 
tance between the two centroids. Note that the orientation of the 
direction connecting the two centroids, in the assumption of a ro- 
tating disk kinematics, identifies the disk line of nodes. Figures 
|2] to [4] show as an example the application of the method to the 
object SSA22a-C16 from the AMAZE sample. 

We can understand the physical meaning of the spectroas- 
trometric characteristic radius (r spec ) in the assumption of rotat- 
ing disk kinematics. We obtain the blue and red images by bin- 
ning the spectral dimension of the data cube, and then measur- 
ing the light centroid position for each of the two resulting bins. 
Therefore we can consider the blue and red centroids, and their 



5 



A.Gnerucci: Dynamical mass estimator based on spectroastrometry. 



Blue image 



Red image 




-2-10 1 2 

Arcsec (referred to the abject position) 




-3-2-10 1 2 3 

Arcsec (referred to the abject position) 



Fig. 3. Example of the blue and red image for the same object of Fig. El White isophotes and black dashed lines represent the best 
model and axis orientation for the two dimensional gaussian function fitted to the images. The blue and red filled circles represent 
the position of the centroid of respectively the other image. 
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Fig. 4. Map of the blue and red centroids for the same example 
of Figs.[2]and[3]The dotted gray grid represents the pixel layout. 



associated average velocities, as a "two-point" rotation curve of 
the gas: each point is defined by its velocity (+cr and -<x) and its 
distance from the center (+r spec and -r spec ), i.e. from the global 
peak of line emission. We note that we are actually measuring 
the rotation curve of the object using spectroastrometry as in 
iGnerucci et al.l (120101) . but with the difference that here we are 
considering only two spectral bins. This is a natural consequence 
of the poor S/N data we have to deal with. 



Finally we can define the spectroastrometric mass estimator 
(M spec ) as: 



M S p ec sin i — fspec ftspec — fs[ 



FWHM 2 r 



spec 



(2) 



where f spec is a calibration factor that takes into account both 
the geometry of the mass distribution and our choices for mea- 
suring velocities and dimensions. The value of f spec can be cali- 
brated empirically by comparing the virial product with the more 
accurate masses from full dynamical modeling. We have inserted 
the factor sin 2 i in the first term because of the coupling to the 
disk inclination discussed above. 

The principal advantage of our estimator with respect to the 
classical virial mass is the way we estimate the characteristic di- 
mension of the galaxy. Firstly, because we estimate it directly 
using the ionized gas emission, it is more directly related to the 
observed kinematics than the continuum. In addition, its mea- 
surement is not affected by the finite spatial resolution, and we 
can accurately estimate sizes down to ~ 1/10 of the spatial res- 
olution. 

As previously observed, this is the principal advantage of the 
spectroastrometric method and becomes a fundamental feature 
when dealing with objects observed with poor intrinsic spatial 
resolution as in the case of high z galaxies. In the example pre- 
sented in Figs. [2] [3] and [4] we can estimate the blue and red 
centroid positions, and consequently r spec , with an accuracy of 
~ 0.03" that is ~ 1/4 of the pixel size (0.125") and ~ 1/20 of 
the spatial resolution of the data (for the data used in the example 
and for most of the data used in this paper the spatial resolution 
is ~ 0.6" FWHM). With any other method, the estimate of the 
galaxy half light radius needs to be corrected for beam smearing 
with consequently high uncertainties especially when the object 
is only marginally resolved. 

Another important point is that cr and r spec encode different 
pieces of information about the gas kinematics but are extracted 
from the same data, thus avoiding possible biases which might 
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affect the dynamical mass estimate when, e.g., using the contin- 
uum emission to estimate sizes. We conclude this section with a 
remark: as for the classical virial mass estimator and dynamical 
mass estimates from full dynamical modeling, with our spec- 
troastrometric virial mass estimator we implicitly assume that 
the gas kinematics is that of a rotating disk. As such, we can ap- 
ply it only to objects that show gas kinematics consistent with 
rotation, or which are spatially unresolved. 

5.1. Calibration of the spectroastrometric estimator 

The next step is to calibrate the spectroastrometric mass estima- 
tor and verify its reliability with real data. We thus compare the 
estimator predictions with the more accurate estimates from full 
dynamical modeling of ourz ~ 3 galaxies from the AMAZE and 
LSD samples and of the z ~ 2 subsample of h igh S/N rotating 
disks from the SINS survey dCresci et al.ll2009l) . 

In the right panel of Fig. [TJ we compare the spectroastromet- 
ric virial product (p. spec - FWHM 2 r spec /G) and the dynamical 
mass (Mdyn) resultin g from complete dynam ical modeling of the 
galaxies analyzed in iGnerucci et al.l (1201 ll) . The values for the 
dynamical mass of the best fit models and for the spectroastro- 
metric mass estimator of the various objects analyzed in this pa- 
per are also reported in table [TJ We also selected 8 objects with 
high S/N, "rotating disk" gas kin ematics and robus t dynamical 
modeling from the SINS sample dCresci et al . 2009). There is a 
tight linear correlation between fi spec and Md yn sin 2 i. 

By fitting the linear relation to the data 

\og{M dyn sin 2 i) = logQi spec ) + log(f spec ) (3) 

we obtain log f spec = 0.02 + 0.05,1.6./,^ = 1.0*!H with a resid- 
ual dispersion of 0.15 dex. The right panel of Fig. [TJ should be 
compared with the left panel where we have plotted the simi- 
lar relation but using the "classical" virial product //,,,> instead of 
the spectroastrometric virial product fi spec . With our new dynam- 
ical mass estimator the relation with M dyn sin 2 i becomes tighter 
(the rms of the residuals is decreased by 0.15 dex) and, most 
importantly, becomes linear compared to the non-linear relation 
with the classical virial estimator. In conclusion, compared to 
spectroastrometric estimates, "classical" virial masses are more 
biased by systematic effects (especially in the low mass range) 
and are less accurate proxies of the dynamical mass. 

Finally, using the fitted calibration factor, we can obtain the 
value of the inclination-independent dynamical mass for an ob- 
ject based on the spectroastrometric mass estimator as 

M dyn sinH = l.O + _o\n spec (4) 

with a systematic error of +0.15 dex. 

We have performed a set of simulations in order to under- 
stand how the spectroastrometric mass estimator is affected by 
various parameters (i.e. the object's own dynamical features or 
the instrumental setup) and to confirm the results obtained with 
our data sample. These simulations, which are presented in de- 
tail in AppendixfA] indicate that 

- the spectroastrometric dynamical mass estimator is indepen- 
dent of the spatial resolution of the data (within the "observa- 
tional" errors) for a variation of the PSF FWHM in the range 
0.1" -1.0"); 

- the estimator depends on the disk inclination like the dynam- 
ical mass, following a sin 2 i law; 

- the estimator allows one to recover the dynamical mass with 
a rms dispersion of 0.3 dex when using for the simulations a 



S/N, disk scale length and dynamical mass values similar to 
those of the real data. 
- for an exponential disk model r spec correspond on average to 
a fraction of ~ 0.75 of the disk scale length, although this 
relation depends on the assumed disk profile. 

6. Application of the spectroastrometric estimator 
to the AMAZE and LSD complete samples 

In this section we apply our calibrated estimator to the galaxies 
from the AMAZE and LSD samples where a dynamical mass 
estimate from full dynamical models has not been possible. Our 
estimator can be applied only to objects with rotating disk kine- 
matics. Therefore we exclude all galaxies which are spatially 
resolved and show complex kinematics inconsistent with a rotat- 
ing disk. Clearly, we include all the objects which are unresolved 
and/or possess low S/N because in those cases we have no hints 
about their kinematics. This selection is the same that would be 
applied if using the classical v irial method. 

As previously observed in IGnerucci et al.l (1201 ll) we imple- 
mented a simple but quantitative method to identify objects with 
gas kinematics consistent with disk rotation based on classifying 
a galaxy as "rotating" or "non-rotating" according to whether its 
velocity map is well fitted by a plane. Using this method we iden- 
tified a subsample of 1 1 "rotating" objects and 6 "unclassified" 
objects. Therefore we apply the spectroastrometric estimator to 
these "rotating" and "unclassified" objects. 

We first calculate fi spec and then convert it to M dyn sin 2 i as 
explained in the previous section. We assume an average disk 
inclination of i = 60° (the mean inclination over a population of 
randomly oriented disks) to finally obtain the dynamical mass. It 
was not possible to derive the dynamical mass estimator M spec 
for all objects. In same cases the S/N of the spectrum was too 
low, in other cases the galaxies were extended and too clumpy 
for a reliable estimate of the red and blue centroids. Finally, 
we obtain a total of 16 dynamical mass estimates (14 from the 
AMAZE sample and 2 from the LSD sample) ranging from from 
2.5 x 1O 7 M to 2.8 x 10 10 M o with a mean value of 1 .3 x 1O 9 M . 

The values obtained are presented in figure [5] where we 
compare our spectroastrometric mass estimates with the stel- 
lar masses of the "rotating" or "unclassified" objects in the 
AMAZE and LSD samples (we use the spectroastrometric esti- 
mates even for the galaxies where a full dynamical estimate was 
possible). The stellar masses for this subsample of the AMAZE 
and LSD samples are es timated from standa rd broad-band SED 
fitting (see Sect. 6 o f iMaiolino^ t al. 200 8ai and Sec t. 4.3 of 
iMannucci etai1l2009l) . and scaled to a lchabrierl d2003l) IMF di- 
viding flieM i _values calculated using a Salpeter IMF for a factor 
~ 1.7 dPozzetti et al.ll2007l) . We note that we estimate M dyn sin 2 i 
and assume i = 60deg. Therefore, in comparing spectroastro- 
metric and stellar masses we have to take into account the effect 
of the disk inclination on the dynamical masses. In figure [5] we 
plot vertical bars for the spectroastrometric dynamical mass cor- 
responding to inclinations of 90°, 60°, 10° (lower, central and 
upper value respectively). 

The comparison between dynamical masses for i = 60° and 
stellar masses, indicates that about one fourth of the dynamical 
masses are unphysically underestimated with respect to stellar 
masses. This fraction is mainly composed of "unclassified" ob- 
jects. Indeed, 2 out of 5 "unclassified" objects have M d y n lower 
than M, for more than 0.2dex while this is true for only 2 out of 
the 1 1 "rotating" objects. 

Apart for the non-negligible systematic errors affecting stel- 
lar masses, the causes behind an underestimate of dynamical 
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UDjecr ' 


Sample 




z 


r WMM ykms J 




R c ["] (3) 




(5) 

r 


^spec 


r"i(4) 


(5) 

h*spec 


fl/f n.'«2 J (5) 

M j vtt sin i y ' 


SSA22a-M38 


AMAZE 


3 


288 


489 ± 36 





J 1 


± U.UJ 


11.38 


±0.06 


0.29 


±0.02 


11.08 ±0.16 


11.33 ±0.02 


SSA22a-C16 


AMAZE 


3 


065 


274 ± 24 





3 J 


± u.uz 


10.82 


±0.08 


0.09 


±0.01 


10.05 ± 0.20 


10.19 ±0.02 


CDFS-2528 


AMAZE 


3 


689 


295 ± 47 





1. 1 
J 1 


± U.Uo 


10.83 


±0.13 


0.13 


±0.02 


10.24 ± 0.3 


10.59 ± 0.05 


SSA22a-D17 


AMAZE 


3 


079 


514 ± 108 





17 


±0.11 


11.24 


±0.28 


0.08 


±0.01 


10.57 ±0.43 


10.52 ± 0.06 


CDFA-C9 


AMAZE 


3 


219 


320 ± 25 





19 


±0.09 


10.85 


±0.22 


0.042 


± 0.003 


9.87 ±0.16 


9.81 ±0.03 


CDFS-9313 


AMAZE 


3 


652 


264 ± 25 





22 


±0.09 


10.68 


±0.19 


0.016 


± 0.009 


9.27 ± 0.58 


9.14 ±0.04 


CDFS-9340 


AMAZE 


3 


652 


150 ±46 





19 


±0.12 


10.18 


±0.30 


0.09 


±0.03 


9.55 ±0.56 


10.1 ±0.1 


3C324-C3 


AMAZE 


3 


282 


303 + 51 





50 


±0.06 


10.99 


±0.09 


0.09 


±0.02 


10.13 ±0.42 


9.84 ±0.12 


CDFS-14111 


AMAZE 


3 


596 


217 + 26 





21 


±0.08 


10.47 


±0.17 


0.046 


± 0.007 


9.56 ± 0.24 


9.59 ± 0.04 


CDFS- 16767 


AMAZE 


3 


615 


281 ±44 





16 


±0.10 


10.63 


±0.27 


0.05 


±0.01 


9.80 ±0.34 


9.76 ±0.11 


Q0302-C131 


LSD 


3 


240 


176 ± 22 





07 


±0.02 


9.93 


+ 0.38 


0.044 


± 0.007 


9.36 ± 0.24 


8.98 ± 0.07 



Notes. The spatial resolution of observations was ~ 0.6" (FWHM). (2) Corrected for the instrumental spectral broadening by subtracting in 
quadrature the spectral resolution. (3) Corrected for the instrumental beam smearing by subtracting in quadrature the PSF FWHM. {4) No correction 
for instrumental beam smearing, the spectroastrometry technique enables to make distances measurements down to ~ 1/ 10 of the spatial resolution 
(see text). (5) Masses are in unity of log(M/Mo). 



masses should be ascribed either to the unknown disk inclination 
and/or to the possible contribution of non-rotational motions in 
"unclassified" objects. 

Regarding the first effect, we assume a mean inclination of 
60° in our estimates, but we cannot verify this assumption. The 
vertical bars plotted in the figure identify the 10 - 90deg in- 
clination range and show that, when taking into account incli- 
nation effects, only one object is not consistent with the condi- 
tion M* < Mdyn- Regarding the second effect, we note that the 
largest discrepancies are for the "unclassified" objects for which 
the assumption of rotating gas disks cannot be verified. Indeed 
the presence of non-rotational motions could significantly affect 
the spectroastrometric mass estimate. We will discuss in more 
detail all the effects that can bias our dynamical mass estimate 
in the next section. 



7. Discussion and conclusions 

We have presented a new method to estimate galaxy dynamical 
masses based on the spectroastrometry technique which is an 
improvement over "classical" virial mass estimates. This method 
is potentially very useful for its application to high z galaxies 
where the intrinsic spatial resolution and signal to noise ratio of 
the observations is often too poor to allow for the more reliable 
analysis based on spatially resolved kinematical modeling. 

The principal advantage of spectroastrometry is its capability 
of overcoming the limit of the spatial resolution of the observa- 
tions and this feature is important in our application of the tech- 
nique. We estimated the characteristic dimensions of line emis- 
sion regions down to 0.04" with a > 3cr significance; since the 
average spatial resolution is ~ 0.6", the spectroastrometry has 
allowed us to directly probe spatial scales down to ~ 1/15 of 
the spatial resolution. Our measurements do not rely on the de- 
convolution techniques which are used in the more conventional 
analysis. 

We have shown that our spectroastrometric dynamical mass 
estimates correlate well with dynamical mass estimates from full 
dynamical modeling, better than classical virial mass estimates 
(Figs. [2] and 0. We can conclude that the characteristic radius 
inferred with spectroastrometry has not only a better accuracy, 



12 r 




8 9 10 11 12 

log lo (M./M ) 



Fig. 5. Spectroastrometric vs stellar mass for AMAZE and LSD 
objects. Filled circles and open diamonds corresponds to ob- 
jects classed as "rotating" and "unclassified" respectively, and 
the spectroastrometric masses have been computed for a disk in- 
clination of 60 deg. The vertical bars show the effect on M spec 
of extreme disk inclinations of 90° and 10° (lower and upper 
values respectively). The gray area is the region of unphysical 
M+ > M spec . Brown dashed lines represents respectively the 
30%, 60% and 90% gas fraction loci. Only the "unclassified" 
objects have dynamical masses {M spec ) significantly lower than 
the corresponding stellar masses. 

but it is also more reliable compared to the "classical" estimate. 
An additional advantage is that a continuum image of the galaxy 
is not needed, but it is possible to measure velocities and dimen- 
sions only from gas emission lines. 

There is a fundamental requirement/assumption for spec- 
troastrometric dynamical mass estimates as well as for classical 
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virial ones: the gas kinematics must be that of a rotating disk. 
Only then is the shift between the red and blue image light cen- 
troids connected with disk rotation. In fact, if the gas kinematics 
is that of a rotating disk, the blue- and the red-shifted gas emis- 
sion come from spatially distinct regions (i.e. the two halves of 
the disk) and the corresponding light centroids will be shifted. 
In contrast, for a non-rotating disk, the blue- and red-shifted gas 
emission could be spatially mixed, the shift between the respec- 
tive centroids will be smaller and the connection with the galaxy 
dynamical mass via the rotating disk model would no longer be 
possible. Of course, the assumption of a rotating disk needs to 
be valid also to have a connection between the observed line 
width and the galaxy dynamical mass. In this case, the line width 
is simply due to unresolved rotation. Conversely, any important 
non-rotational contribution to the line width will spuriously in- 
crease our mass estimate. This effect is present and unavoidable 
also in "classical" virial mass estimates. 

To obtain a reliable dynamical mass estimate, it is important 
to verify that the object exhibits rotating disk gas kinematics. 
For the pur poses of this paper we adopted the simple method 
presented in lGnerucci et alj d201 lb based on the^ 2 statistics to 
identify smooth velocity gradients in the line-of-sight velocity 
map of the object (a smooth velocity gradient is in fact, for this 
low S/N and poor spatial resolution object, a good signature of a 
rotating disk kinematics). We remark that this simple method is 
used principally to discard objects inconsistent with disk rotation 
rather than to analyze the consistent ones. In fact our spectroas- 
trometric mass estimator is mainly useful for low S/N and poorly 
resolved objects and in such cases even a good measure of the 
X 2 statistic to assess the presence of a smooth velocity gradient 
is difficult to obtain. 

The easiest application of the method is for objects which do 
not show a complex morphology (e.g. with many knots of emis- 
sion). This is often the case for high-z galaxies which are mostly 
barely resolved. When the object is well resolved and the kine- 
matics is that of a rotating disk, then a full modeling is possible. 
Therefore the spectroastrometric mass estimator is ideal to apply 
to sources which, for lack of S/N or spatial resolution, cannot be 
analyzed with standard kinematical models. 

To calibrate the spectroastrometric mass estimator and verify 
its reliability we applied it to 19 high redshift galaxies, whose 
SINFONI spectra had been already analyzed and modeled: 1 1 
Z ~ 3 objects from the AMAZE and LSD samples and 8 z ~ 2 
objects from the SINS sample. For each galaxy we compared 
the estimator prediction with the dynamical mass value resulting 
from the spatially resolved kinematical modeling. We observed 
a tight correlation between our estimator prediction and model 
dynamical mass values (dispersion of 0.15 dex), indicating that 
the spectroastrometric estimator is reliable. From this compari- 
son we also calibrated it, obtaining 

(±0.15 dex - systematic) (5) 

where V is the line FWHM, and r spec is half the distance be- 
tween the red and blue light centroids. Such calibration, and the 
confidence of the reliability of the estimator, can of course be 
strengthened by enlarging the sample of galaxies for which both 
the spectroastrometric and standard dynamical mass are mea- 
sured, extending the range of mass and redshift of the sources. 
The classical virial mass estimates are a much worse proxy of 
the true dynamical mass (non-linear relation and 0.3 dex disper- 
sion). We remark that we used 3 different samples of objects at 



different redshifts for this calibration, obtaining a tight correla- 
tion over a full two orders of magnitude in Md v „sin 2 i. 

As an example of its application, we finally applied our es- 
timator to 16 objects of the AMAZE and LSD samples show- 
ing that it is possible to estimate a wide range of masses from 
1 .3 x 1 8 M to 1 . 5 x 1 1 1 M sm with a mean value of 7 . 1 x 1 Q 9 M mn . 
We then compared spectroastrometric and stellar masses, Md y „ 
vs. M*, showing the derived masses are consistent with the con- 
dition M* < Md v „, once the (unknown) disk inclination is taken 
into account. Indeed, the estimator presented in the above equa- 
tion is actually M^sin^i , and the inclination dependence can 
be corrected by assuming an average value of i = 60 deg, if no 
other information is available. 

Our spectroastrometric estimator is given by the product of 
the spectroastrometric radius r spec with the line FWHM squared, 
therefore the principal biases to our mass estimate will lie in 
these two quantities. 

A first bias can originate from the unverified assumption of 
rotating disks. In this case, as previously observed, the r spec es- 
timate will be smaller because of the spatial mixing of blue and 
red-shifted gas emission leading to a underestimated dynamical 
mass. In our z ~ 3 sample, the low S/N of data can be responsi- 
ble, in some cases, for the selection of intrinsically non-rotating 
galaxies (especially for the "unclassified" objects) introducing 
spurious low-mass points in the plot of Fig. [3] 

Another effect that could bias the dynamical mass estimate is 
the presence of low-inclination (face-on) rotating disks. In this 
case the line FWHM will be lower because of the projection 
along the line of sight. Also r spec could be biased toward small 
values because the lower surface brightness of the outer regions 
of a face-on disk can push the estimated position of the blue and 
red centroids toward the disk center. Also in this case we will 
obtain a lower dynamical mass estimate mainly because of the 
assumed inclination of 60°. We note that in a high-z sample like 
ours, the cases with lower S/N and spatial resolution (i.e. "un- 
classified" objects) can correspond to low inclination objects. In 
fact lower surface brightness and a velocity gradient consistent 
with zero can be indications of a face-on disk. Hence an inclina- 
tion smaller than 60° might be more representative for these ob- 
jects. Unfortunately it is very difficult to assess this effect more 
quantitatively. 

Nonetheless it is important to observe that in the case of a 
low mass estimate it is not possible to distinguish between the 
above mentioned effects because of the previously noted kine- 
matical ambiguity for "unclassified" objects (i.e. we cannot dis- 
tinguish between low-inclination disks or non-rotating object). 

Another effect that can bias our dynamical mass estimate is 
the presence of an important contribution to the line width from 
non-rotational motions. This effect rises the mass estimate and 
can also produce minor biases on the r spec estimate. In fact, if the 
line profile is broadened by non rotational motions, we would 
insert spurious high velocity planes from the line wings in the 
"red" and "blue" images of the source that can alter the centroid 
position in an unpredictable way. Many authors suggests that 
the observed high turbulence of high z disks can provide further 
dynamical support and add to the dynamical mass a contributio n 
of "sigma-supported" mass (see for example fepinat et al1 l2009). 
However, in our method it is not possible to identify such cases 
and quantify this effect because of the S/N and spatial resolution 
of our data and because the principal effect of a broadening of 
the line from non-rotational motions is to raise the dynamical 
mass estimate. This effect is present and undetectable also in 
"classical" virial mass estimates. 
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Appendix A: Simulations 

We present here a set of numerical simulations in order to un- 
derstand how the spectroastrometric mass estimator is affected 
by several parameters such as a galaxy's dynamical structure or 
the instrumental setup. 

Our simulations are based on a dynamical model of a high z 
galaxy. The model consists of a rotating thin gas disk for which 
we neglect hydrodynamical effects, therefore the disk motion 
traces the full gravitational potential. The gravitational poten- 
tial is given by an exponentially-decaying disk-like mass dis- 
tribution with scale radius ro and a total mass Md y „ at a lOkpc 
radius (this i s the same definition of dynamical mass given in 
ICresci et~al1 (|2009)). The galaxy is placed at redshift z and has 
an inclination i with respect to the line of sight. We simulate in- 
tegral field spectroscopic observations of a gas emission line by 
setting the principal instrumental parameter (spatial and spectral 
resolution and x, y and wavelength binsize) to match those of 
SINFONI, the instrument with which the data presented in this 
paper has been obtained. F or a detailed descriptio n of the model 
refer, e.g., to appendix B of M arconi et alj d2003l) . 

We simulate the presence of noise by adding to the resulting 
data cube normally distributed random numbers, characterized 
by zero mean and standard deviation <x„ 0(> . This is chosen such 
that spectra of the galaxy, extracted from apertures similar to 
those used with real data, have the same S fN as observations. 

Finally, we analyze these simulated data cubes in the same 
way as the real ones and calculate the spectroastrometric virial 
product {p spec = FWHM 2 r spec /G). The uncertainties of the 
"simulated" [i spe c values depend on the uncertainties on the mea- 
sured photocenters positions and line FWHM, just as in the real 
data. 

The parameters of the baseline model are chosen to match 
the average characteristics of the observed objects: 

- redshift z = 3, 

- characteristic radius ro = 2kpc, 

- dynamical mass Md yn — 10 n MQ, 

- inclination i = 60°, 

- spatial resolution 0.65"FWHM, 

- spatial scale 0.125" x 0.125", 

- spectral resolution 80£hjs~'FWHM, 

- dispersion 15kms~ l . 

In the following, we will vary some of these parameters and ob- 
serve how this affects the estimated dynamical mass M spec with 
respect to the "true" model dynamical mass Mdyn- The ranges of 
variation of these parameters are chosen to approximately match 
the distribution of the various parameters values observed in our 
galaxy samples. 

In the first set of simulations we test the effect of vary- 
ing the spatial resolution of the observations. In Fig. I A. 1 1 we 
show the ratio (in log scale) between the model dynamical mass 
M c i vn sin 2 i and the spectroastrometric virial product (ju spec - 
FWHM 2 r spec /G) for various spatial resolution values. The val- 
ue s are perfectly consist ent, confirming what we already found 
in iGnerucci et al.l yOlO) that the spectroastrometric technique 
is independent of the data spatial resolution. For the baseline 
model only and the simulation presented in the figure, the ra- 
tio fi spec /Mdy„sin 2 i has a mean value of 0.07dex with a dis- 
persion of the points of ~ 0.03dex. We remark that the cal- 
ibration factor determined in the paper represents the average 
HspeclMd yn sin 2 i for the analyzed galaxies. 

In the second set of simulations we test the effect of varying 
the disk inclination, but keeping all the other parameters of the 
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Fig.A.l. Simulations. Ratio between the spectroastrometric 
virial product fi spec = FWHM 2 r spec /G and the model Md y „sin 2 i 
for simulations with various PSF FWHM values. The dashed line 
represents the mean value of the ratio. Dotted lines represent the 
rms dispersion of the points around the mean value. 



basic model fixed. We expect that the spectroastrometric virial 
product will vary following a sin 2 i law due to the fact we are 
measuring the velocity of the gas along the line of sight; there- 
fore, the ratio fi spec /Md v „sin 2 i should be independent of incli- 
nation. In Fig. I A. 2 1 we show the dependence of this ratio with 
the inclination angle. p sp< , c /Md y „sin 2 i is systematically under- 
estimated by ~ 0.2 - 0.4dex in the simulations with i < 45°. 
This discrepancy is due to the fact that the line FWHM be- 
comes smaller than the instrumental spectral resolution at low 
inclinations and the line profile is spectrally unresolved, thus 
affecting the spectroa strometric measurements. As discussed in 
IGnerucci et al.1 d2010l) . the line profiles should always be well 
spectrally resolved to apply spectroastrometry. Therefore, we 
repeated the simulations at low-inclination increasing by a fac- 
tor 2 both the spectral resolution and dispersion with respect to 
the baseline model: the new jj spec /Md yn sin 2 i values are marked 
by the upward arrows in the figure. Once all the simulations 
have been computed with the proper spectral resolution and dis- 
persion, the jj spec /Md Y „sin 2 i ratio remains fairly constant as ex- 
pected. We finally remark that the problem of poor spectral reso- 
lution is not present in our data, since the line FWHM are always 
much larger than the instrumental spectral resolution. 

In the third set of simulations we vary the disk scale radius ro 
and the dynamical mass Md y „. ro varies in the range (0.5 - 8 kpc) 
and Mdyn m the range (10 9 - 10 12 Mq)- The values for the simu- 
lations have been selected following a comparison with the ob- 
served (ro, Mdyn) distribution for our sample, which is shown in 
Fig. lA.3l (we plot the observed (ro, Md yn ) values as open triangles 
and the values selected for simulations as filled circles). We note 
that ro is the ionized gas scale length and Md yn is the dynamical 
m ass from the ful l dynam ical modeling of the objects analyzed 
by IGnerucci et al.l d201 lb and used in this paper. 



In Fig. IA.4l we show the ratio fx spec /Mdy n sifi i as a function of 
Mdyn (left panel) and ro (right panel). In each panel we vary the 
color of the symbols to identify either ro (left panel) or Md yn 
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Fig. A.4. Simulations. Left panel: measured fj. spec /Md yn sin 2 i values from simulations obtained with various combinations of Md yn 
and ro, plotted as a function of Mdyn- Colors denote different ro values. Right panel: same as left panel but fi spec /Md yn sin 2 i is now 
plotted as a function of ro and colors denote different M</,.„ values. In both panels gray triangles represent the measurements from 
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Fig.A.2. Simulations. Ratio fi spec /(Md v „sin 2 i) for simulations 
with various disk inclination values. At i < 45° the line profiles 
become unresolved, thus violating the conditions under which 
spectroastrometry can be applied. To remedy this, the upward 
arrows at i = 20° and 30° indicate measurements for which the 
spectral resolution and dispersion were increased. The dashed 
line represents the mean value of the ratio. Dotted lines repre- 
sent the rms dispersion of the points around the mean value. 



(right panel). We also plot for comparison the measurements 
from the real data (gray triangles). Figure clearly shows that our 
simulations well reproduce the real measurements with a mean 
ratio HspeclMjvnSin 2 ! of 0.05dex and a dispersion of ~ 0.3dex 
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around this value, consistently with the f spec value obtained in 
Sec. O 

Finally we compare the measured spectro astro me trie char- 
acteristic radius r spec with the model disk scale radius ro. In 
Fig. IA.5l we show this comparison for all the models plotted in 



Figs. IA.3l and lA.4l From figure we can observe that r spec is sim- 
ilar to the disk scale radius ro for the smaller disks (ro < 2kpc) 
whereas for larger disks it is smaller than ro by ~ 20% - 30%. 
From the definition of fi spec (eq. |2]) we expect that this ap- 
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parent underestimate of the disk scale radius can affect in the 
same sense the fj. spec proxy of the dynamical mass. Indeed from 
Fig. IA. 41 we can not observe any evident negative bias for larger 
disk scale radii. This is a clear indication that the size needed in 
combination with the line width to obtain the mass is not accu- 
rately given by ro, but by our spectroastrometric radius. 
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Fig. A.5. Simulations. Measured r spec values for all the models 
of Figs. IA.3l and lA.4l as a function of model disk scale radius ro. 



Finally, based on the simulations presented we can conclude 



that 



- the spectroastrometric dynamical mass estimator is indepen- 
dent of the spatial resolution of the data (within the "observa- 
tional" errors) for a variation of the PSF FWHM in the range 
0.1" -1.0"); 

- the estimator depends on the disk inclination like the dynam- 
ical mass, following a sin 2 i law; 

- the estimator allows one to recover the dynamical mass with 
a rms dispersion of 0.3 dex when using for the simulations a 
S/N and (ro, Md v „) values similar to those of the real data. 

- for an exponential disk model r spec correspond on average to 
~ 0.75ro, although this relation depends on the assumed disk 
profile. 
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